Determining Pair Interactions from Structural Correlations 
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We examine metastable configurations of a two-dimensional system of interacting particles on a 
quenched random potential landscape and ask how the configurational pair correlation function is 
related to the particle interactions and the statistical properties of the potential landscape. Un- 
derstanding this relation facilitates quantitative studies of magnetic flux line interactions in type II 
superconductors, using structural information available from Lorentz microscope images or Bitter 
decorations. Previous work by some of us supported the conjecture that the relationship between 
pair correlations and interactions in pinned flux line ensembles is analogous to the corresponding 
relationship in the theory of simple liquids. The present paper aims at a more thorough under- 
standing of this relation. We report the results of numerical simulations and present a theory for 
the low density behavior of the pair correlation function which agrees well with our simulations and 
captures features observed in experiments. In particular, we find that the resulting description goes 
beyond the conjectured classical liquid type relation and we remark on the differences. 



I. INTRODUCTION 

Magnetic flux lines in type II superconductors be- 
long to a class of systems in which interacting particles 
are subjected to random pinning forces @,§. At very 
low temperatures, the resulting configurations are deter- 
mined by a competition between elastic energy cost and 
pinning energy gain as the flux line ensemble attempts 
to accommodate the random medium. This competi- 
tion gives rise to a rich and complex energy landscape 
of metastable states. One expects that structural prop- 
erties of the particle configurations should contain infor- 
mation about the interparticle interactions as well as the 
statistical properties of the pinning potential landscape, 
the pinscape. 

In an earlier publication || by some of us, henceforth 
referred to as SHTCG, the near-equilibrium configura- 
tions of flux lines in a Nb thin film obtained from Lorentz 
microscope images [Q were analyzed. SHTCG argued 
that a sequence of flux line configurations on the fixed 
pinscape resembles the instantaneous configurations of a 
classical simple liquid, with the underlying pinscape play- 
ing the role of an effective heat bath. This conjecture 
implies that in the low-density limit, the pair-correlation 
function g(r) should be of the form || 

g(r) = e-W- + . . . , (1) 

where U (r) is the flux line pair interaction, U a is an ef- 
fective pinning energy scale, and the terms omitted cor- 
respond to higher order corrections due to many-body 
effects. For stiff magnetic flux lines, the interparticle po- 
tential U(r) (per unit length) is given in the low density 
limit by the asymptotic form |].|7j 



where <p Q is the flux quantum and A is the London pen- 
etration depth. By treating U as an adjustable param- 
eter, SHTCG showed that the small r behavior (lift-off) 
of the pair correlation function agrees well with the form 
given in Eq. ([!]), yielding a value for the penetration 
depth in good agreement with the accepted value for Nb 
(see Fig. |l|). Bitter decoration images of flux line configu- 




U(r) 



(2) 



FIG. 1. Pair correlation function obtained from averaging 
over Lorentz microscope images of flux lines in Nb j^,^] at 
T — 5 K. The region of interest is where g(r) rises from zero 
(lift- off). The inset shows a comparison of the lift-off region 
of g(r) with liquid theory behavior, Eq. (|l]) (refer to section 
II for further details). 



rations in the high temperature superconductor BSCCO 
analyzed in the same manner also yield results in good 
quantitative agreement with the accepted value of the 
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penetration depth for BSCCO gj. The fit value for U a 
presumably describes the characteristic energy scale for 
the pinscape disorder. 

The present paper aims at a more thorough under- 
standing of how structural correlations arise from micro- 
scopic properties such as the particle interactions and 
potential landscape. Particular questions we address in- 
clude: Under what conditions can the static configura- 
tions of an interacting system in a quenched random 
potential give rise to liquid-like correlations of the form 
Eq. ([!])? How does Eq. (|l]) compare to a theory where 
the quenched nature of randomness is accounted for ex- 
plicitly? 

We report the results of numerical simulations and 
present a theory that captures the essential features ob- 
served both in experiments and in our simulations, pro- 
viding thus some answers to the questions raised above. 
In particular our simulations show that a relation of the 
form Eq. (Q) emerges in the limit of point-like pins, i.e. 
pins of small spatial extent, and for broad {e.g. ex- 
ponential) pinning strength distributions. In contrast, 
pinscapes composed of identically strong or spatially ex- 
tended pins are not well described by Eq. ([!]) . Our theory 
is based on the observation that the behavior of the pair 
correlation function at small r is dominated by strongly 
pinned flux lines. In the low-density limit, the functional 
form of the pair correlation function therefore arises es- 
sentially from a convolution of the tails of the pinning 
energy distribution with the probability distribution for 
locating a flux line at a given position. 

The paper is organized as follows. We present our 
model and its motivation in Section II. Section III con- 
tains details of the numerical simulations. Numerical re- 
sults and their discussion are given in Section IV, while 
Section V is devoted to a description of the underlying 
theoretical picture. We conclude with a discussion and 
summary of our results in Section VI. 



II. THE MODEL 

We consider a two-dimensional system of N particles 
in a quenched random pinscape consisting of discrete pin- 
ning sites positioned at random. The particles, which can 
be interpreted as stiff magnetic flux lines, have positions 
that are specified by (two-dimensional) vectors r^. As a 
convention, particle and pinning site locations are sub- 
scripted with Latin and Greek indices, respectively. The 
force on the i th flux line at position is given by 



N 



.)• (3) 



Here the first and second sum are the contributions to the 
total force on a given particle due to its interactions with 
the other particles and with the pinning sites; U(r) is the 



interparticle interaction, whereas V a (r) is the particle- 
pin interaction. For stiff magnetic flux lines, the inter- 
particle potential U(r) (per unit length) is given in the 
low-density limit by the asymptotic form Eq. (||). We 
assume that the pinning potential due to a single pin is 
attractive and short ranged, 



V a {r) = -V a e 



2 /r 2 



(4) 



where V a is the pin's strength, and r p its range. 

The locations of the pinning sites r Q are assumed to 
be distributed at random. We seek stable and static flux 
line configurations {r*} given by solutions to 



Fi(rJ,r|,. 



0, 



(5) 



with the additional stability constraint that the matrix 
— dFi/dr* be positive-definite. If {r*, r%, ■ • ■ , r^} is such 
a configuration, the pair correlation function g*(r) for 
this configuration is given by 



>(r) = 



qy(x-rv(x)dx) 

[fp*(x)dxP 
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*(x-r*), 



(6) 



(7) 



and (. . .)e denoting an average over angles. The pair 
correlation function g(r) is obtained by averaging Eq. (|^) 
over independent configurations. 

If the analogy to liquid structure theory adequately de- 
scribes the emergence of structure in this quenched ran- 
dom system, then Eqs. (|l|) and imply that (to lowest 
order) 



ln(— r 2 lng(r)) oc — — . 

A 



(8) 



Thus liquid theory predicts that a plot of the correspond- 
ing quantity versus r should yield a linear relation with 
slope 1/A, as observed experimentally in the inset to 
Fig. | 

Our numerical work involves solving Eq. (|s|) for static 
configurations and thereby obtaining the configuration- 
averaged pair correlation function. Details of the imple- 
mentation and presentation of our numerical results are 
given in the following two sections. 



III. NUMERICAL IMPLEMENTATION 

In this section we describe our numerical implemen- 
tation. Simulations were carried out on a square region 
with 7250 randomly placed pins. Having constructed a 
random pinscape, we placed N — 65 flux lines at random 
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and solved numerically for stable equilibrium configura- 
tions of Eq. (||) using a three stage procedure: (i) Starting 
from a random initial configuration, the vortex system is 
allowed to relax dynamically according to 



dt 



F i (r 1 ,r 2 , . . .,ry), 



(9) 



until the magnitudes of all forces on the RHS are less 
than a prescribed target value, (ii) Once this has been 
achieved, the resulting configuration is used as an initial 
guess in a Newton-Raphson algorithm to find a solution 
to Eq. (^). Solutions found this way are not guaran- 
teed to be stable with respect to small perturbations of 
the particle locations. Therefore, (iii), a linear stabil- 
ity analysis is carried out to check for stability. If the 
configuration obtained this way does not turn out to be 
stable, stage (i) is resumed with a reduced target value. 
The whole procedure is repeated until a stable configu- 
ration has been found. The combination of a molecular 
dynamics algorithm with a Newton-Raphson scheme, as 
described above, turns out to be far more computation- 
ally efficient than using molecular dynamics alone |J. 

Lengths were chosen so that the square corresponds to 
a 5 x 5 /im 2 field of view in a magnetic field of roughly 
50 Gauss. With this choice, all lengths will be reported 
in ^m. The London penetration depth was taken to be 
A = 0.04 /im. These choices resemble the experimental 
condition in Ref . || . The parameters varied in the simu- 
lations were the range of the pinning forces r p , with val- 
ues r p € [0.015, 0.05] /im, and the distribution of pinning 



strengths V a , cf. Eq. (§). 
distributions: 



We considered the following 



V{V a ) = 5{V a -V ), 
for identical pinning strengths, and 

v v {v a ) = ^ e - {c ^ r ,V a >0, 



(10) 



(11) 



for stretched-exponentially distributed pinning strengths, 
with the constants ci and C2 given by c\ = vc2/T(l/ v), 
and C2 = r(i/^) ' wnere ^( x ) ^ s f ne Gamma function. 
In particular, the distributions used were v = 2 (half- 
gaussian), v — 1 (exponential), and v = 1/2 (stretched 
exponential). The mean of the above distributions is 
given by V . The unit of energy was chosen such that 



8tt 2 A 2 V 



200. 



(12) 



This choice is consistent with experimental observations 
H , and ensures that the maximum force exerted by a pin 
is independent of the range of the pinning potential. 

The numerical integration in stage (i) was carried out 
using a 4 th order Runge-Kutta method [|l0| with adap- 
tive stepsize and an initial force resolution of 10~ 6 in 



the above units with subsequent decrements by factors 
of 10, as needed. Algorithms of the software libraries 
LINPACK and EISPACK Q were used in putting to- 
gether stages (ii) and (iii). The maximum residual force 
on any particle in a stable configuration was found to be 
smaller than AF = O(10" 12 ). 

Free boundary conditions were employed, and flux lines 
leaving the square (leaks) were put back in random lo- 
cations. Typically, only a few leaks occurred in a given 
run. 

In order to speed up computations further, several ad- 
ditional measures were taken: If a stable configuration 
was not found after a fixed number of integration steps 
(5000), a new initial configuration was generated. Initial 
configurations were such that flux lines were randomly 
positioned, but no two were allowed to be closer initially 
than 0.2 /im. These measures were found to reduce the 
computation time without changing the results signifi- 
cantly. The particle-pin interaction, Eq. (|]), was cut off 
at a distance of 4r p . Prior to each run, the area was sub- 
divided into a mesh of squares and a look-up table was 
generated providing each square with the location and 
strength of pins within the effective interaction distance 
4r p . This eliminates the need to search all pinning sites 
for each computation of the pinning forces. 
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FIG. 2. Pair correlation function g(r) obtained from aver- 
aging over 212 stable flux line configurations on a pinscape, 
for A = 0.04 and r p = 0.015. The pinning strengths are expo- 
nentially distributed. The inset shows a comparison with the 
liquid theory prediction, Eq. (ph. 



IV. NUMERICAL RESULTS 

In this section we present the results of our simulations 
described above. Fig. || shows the scaling of the pair cor- 
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relation function obtained from averaging over 212 sta- 
ble flux line configurations on a pinscape generated with 
exponentially distributed pinning strengths according to 
Eq. (|n]), with r p = 0.015. The inset to Fig. | shows the 
scaling behavior of ln(— r' ln<7(r)) with r. The solid line 
has slope 1/A. Note the good agreement with the liquid 
theory prediction, Eq. (||) , for values of r in the initial lift- 
off of g(r). The deviations for larger r arise from excess 
correlations (pile-up) near the peak of g(r) not accounted 
for in the naive scaling Ansatz. Such excess correlations 
also are seen experimentally 

Having thus observed numerically a liquid-like rela- 
tion of the form Eq. (||), we next ask how the behav- 
ior of the pair correlation function depends on the range 
of the pinning potential r p . Fig. || shows the scaling of 
ln(— 7*2 lnfif(r)) with r for different values of r p , ranging 
from 0.015 /mi to 0.05 /im. The asymptotic form, Eq. (]§), 
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FIG. 3. Dependence of the pair correlation function on the 
range of the pinning potential r v for exponentially distributed 
pinning strengths. Solid lines correspond to a liquid theory 
scaling, Eq. (^|), with A = 0.04 /im. The values for r p 
are 0.05 /Ltm (a), 0.025 fim (b), 0.017 /Ltm (c), and 0.015 fim 
(d). Averages were calculated from 596, 600, 323 and 212 
configurations, respectively. 



is realized only in the limit of a short-ranged pinning po- 
tential (r p <C A). 

The dependence of the pair correlation function on the 
breadth of the pinning strength distribution is shown in 
Fig. ||. Here we compare the results for configurations 
on a pinscape of spatially narrow pins whose strengths 
were drawn from increasingly broad distributions: delta 
function [Eq. @], u = 2,l and 1/2 [Eq. (§]])]. The pair 
correlation functions corresponding to the three broad 
distributions, v = 2, 1 and 1/2 agree with the predic- 
tion of the liquid-like theory, Eq. (|l|), at least for small r. 
In contrast, Eq. (|l|) does not describe the case of identi- 
cal pinning strengths; the rapid fall-off of ln(— r 2 In g(r)) 
indicates a steeper than expected lift-off of g(r). 
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FIG. 4. Dependence of the pair correlation function on 
the distribution of pinning strengths: (a) identical pinning 
strengths, (b) half-gaussian, v — 2, (c) exponential, v = 1, 
and (d) stretched exponential, v — 1/2. Solid lines indicate 
the scaling, Eq. (^), for A = 0.04 (r p = 0.015). Averages were 
obtained from 148, 228, 212, and 152 configurations, respec- 
tively. 



Thus the results of our numerical simulations clearly 
show a liquid- like behavior of the form Eq. (||) as observed 
in experiments. Furthermore our simulations indicate 
that this scaling only arises in the limit of point-like pins 
drawn from a broad distribution of pinning strengths. 
Its successful application to experimental data suggests 
that the pinscape in these materials also is characterized 
by point-like pins with a broad distribution of strengths. 
Moreover, these observations serve as a warning that sim- 
ulation results based on identical pinning strengths may 
differ subtly from the behavior of physical systems. 



V. THEORETICAL DESCRIPTION 

The analogy drawn between pinned flux lines and ther- 
mally disordered fluids yields an empirically satisfactory 
description of flux line pair correlations. To this point, 
though, it has been justified by heuristic arguments. This 
section introduces a statistical treatment of structure in 
strongly pinned flux lines which provides a superior de- 
scription of disordered flux line distributions yet yields 
the liquid structure result, Eq. (Q), as a limiting case. 

We begin by assuming that the dominant contribution 
to the force on a given flux line is due to its nearest neigh- 
bor. Because of the short-range nature of the flux line 
interactions, this assumption is justified in the limit of 
low flux line concentration, where the average flux line 
separation, as well as fluctuations around it, are much 
larger than the London penetration depth. Furthermore, 
our simulations reveal that more than 99% of the flux 
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lines in a stable configuration are pinned, i.e. within a 
distance of r p /v2 of a pinning site. Thus assuming that 
the flux lines are pinned in pairs is reasonable. In par- 
ticular, this approximation is valid for pairs of flux lines 
pinned at distances smaller than average. These pairs 
determine the behavior of the pair correlation function 
near lift-off, i.e. the region where the function rises from 
zero. They occur only if sufficiently strong pins are avail- 
able. In other words, a pair of nearest-neighbor flux lines 
can be pinned at a separation r only if both pins are 
strong enough to overcome the flux line repulsion force, 
f(r) = -dU{r)/dr. 

From Eq. (Q) we see that a pin a of strength V a can 
exert a maximum force given by 



2 V„ 



e r„ 



(13) 



Let P d {F a ) be the probability that a flux line has come 
to rest on a pin of maximal pinning force F a . The cumu- 
lative probability 



i(F) 



P d {F a )dF a . 



(14) 



is the likelihood that the occupied pin is at least as strong 
as F . If we assume that flux lines are equally likely to 
be pinned anywhere in the sample as long as the corre- 
sponding pinning forces are large enough to overcome the 
flux line interaction, it follows that the pair correlation 
function g(r) has the form 



9(r) = g/(/(r)) 



(15) 



This expression should be valid in the low-density limit, 
where only pairwise-pinned flux lines have to be consid- 
ered and screening effects due to the presence of other 
flux line pairs can be neglected. One can readily show 
that the resulting flux line configurations are stable. 

The distribution of dynamically selected pinning 
strengths, Pd(F a ), does not coincide with the correspond- 
ing distribution of a priori available strengths, P(F a ). 
The latter is readily obtained from the pinning energy 
distribution, V v {V a ), using Eq. (|l3|). Fig. [|] compares the 
distributions of dynamically selected and available max- 
imal pinning forces for exponentially distributed (y = 1) 
pinning strengths, Eq. (11). We see clearly that the par- 
ticles preferentially occupy strong pins. This is consis- 
tent with our expectation that an isolated flux line, or 
one that is effectively isolated, is more likely to be in the 
domain of attraction of a strong pin in its vicinity. 

On the basis of this observation, we hypothesize that 
a given flux line dynamically selects the strongest of n 
pins in its immediate vicinity so that 



Qd(F) = 1 - [1 - Q(F)]' 



(16) 
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FIG. 5. Distribution of available and dynamically selected 
maximal pinning forces. Shown here is the case of exponen- 
tially distributed pinning strengths, v = 1. The dashed line 
corresponds to the distribution of available maximal pinning 
forces as obtained from Eqs. ( |l3| ) and (p"l|). Data points cor- 
respond to the distribution obtained by determining the pin- 
ning strength of the pins associated with every pinned flux 
line. The solid line corresponds to Eq. (pp|), as explained in 
the text. 



P d (F) = nP{F) 



f P{F')dF' 
Jo 



n-l 



(17) 



and thus 



Equations (|l5|) - ( |l7| ) and their successful application to 
our simulations are the central results of this article. 

Fig. | shows that P d (F) is very well approximated by 
n = 2, i.e., 

P d (F) = 2P(F) [ P(F')dF', (18) 
Jo 

implying that its cumulative distribution Q d (F) obeys, 

Q d (F) = l-[l-Q(F)} 2 . (19) 

We find that n = 2 also satisfactorily describes dy- 
namically selected pinning in our simulations with half- 
gaussian and stretched exponentially distributed pinning 
strengths, v = 2 and v = 1/2. 

For exponentially distributed pinning strengths, v — 1, 
we obtain 

(20) 



(21) 



where F a = ^/2feV jr v . Thus, 



Q d (F) = 2e- p / F ° ( 1 - -e 



-F/F a 



5 



and, using Eq. (|l5|), 

g (r) = 4e - 2 'M/ F ° 



1 _ L p -f(r)/F a 
2 



(22) 



Analogous predictions can be obtained for the cases i/ = 2 
and = 1/2. 

Fig. U shows comparisons of the theoretically predicted 
scaling of ln(— lng(r)) with the values obtained from 
our simulations for v = 2, v = 1 and = 1/2. These 
comparisons involve no adjustable parameters once n has 
been fixed (Fig. ||), since the force scale F Q was specified 
for each simulation. The results for all three pinning 
distributions are in very good agreement with the form 
predicted by Eq. ( [l5|) for g[r) near lift-off. 
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FIG. 6. Comparison of the behavior of the pair correla- 
tion function near lift-off against simulation results, liquid 
theor y ty pe scaling (dashed lines) and theoretical prediction, 
Eqs. ( |l^ ) and (|l9|, for half-gaussian (circles), exponential (di- 
amonds), and stretched exponential v — 1/2 (boxes) distri- 
bution of pinning strengths. Note that the theoretical predic- 
tions contain no free parameter once n has been fixed (Fig. |^) , 
while the liquid theory predictions are one-parameter fits to 
the form Eq. (hi) . Curves have been vertically offset for clarity. 



The relative success of the present theory compared 
with the liquid structure Ansatz in predicting the behav- 
ior of g(r) near lift-off emphasizes fundamental differ- 
ences between the two classes of systems. In thermal 
systems such as simple liquids, structural correlations 
arise from the interplay of interaction energies and the 
thermal energy scale. Systems such as pinned flux lines 
arc dominated by quenched disorder and select configu- 
rations in which forces are balanced. Nevertheless, the 
two approaches yield comparable predictions for the flux 
line distribution because of the exponential form of the 
pair potential, Eq. (|^). 

Fig. |](a) reveals that a pinscape generated by identi- 
cal pins differs from the more broadly distributed cases 



in that it gives rise to a distinctively steep lift-off in the 
pair correlation function. This behavior is captured also 
by the present treatment. In the low-density limit, the 
theoretical description, Eq. Jl5|), predicts a step-like in- 
crease in g(r) for the smallest separation at which a pair 
of flux lines can be pinned against their mutual repul- 
sion. Fig. ^(a) shows somewhat smoother behavior. An- 
alyzing the configurations contributing to the lift-off of 
g(r) shows that this smoothing is due to many-body ef- 
fects, i.e. flux lines pinned by the interaction with two 
or more close-by flux lines. Comparing this discrepancy 
to the good agreement with the theoretical prediction in 
the case of broadly distributed pinning strengths leads 
us to conclude that, in these cases, the disorder due to 
variation in pinning strength dominates the contribution 
arising from the spatial disorder of the flux line configu- 
rations. 

A liquid-like form for g(r) can be obtained only if the 
distribution of dynamically selected pins is exponential. 
Eqs. (15) and ( |l6|) show that, in the general case, pref- 
actors and different forms of leading order behavior can 
arise depending on the the tails of the (dynamical) pin- 
ning energy distribution. Our results in Fig. || suggest, 
however, that this dependence is weak and will become 
detectable only for very small values of r. This is consis- 
tent with the success of Eq. (|l|) in describing experimental 
data. We also investigated systems of particles interact- 
ing with power law potentials and found that, when the 
interactions fall off sufficiently fast, the pair correlation 
is again of the form Eq. (|l5|) . 



VI. DISCUSSION AND CONCLUSIONS 

We have shown that the pair correlation function aris- 
ing from stable configurations of interacting particles in a 
random pinscape can contain information about both the 
particle interactions and also the statistical properties of 
the underlying pinscape. Our analysis shows that this is 
the case in the regime where the pinning and pair inter- 
actions are of comparable strength. If the pinning com- 
pletely dominates pair interactions, the pair correlation 
function will essentially reveal the spatial disorder of the 
active pinning sites. In the opposite limit, in which the 
pair interactions dominate, the pair correlation function 
will contain features corresponding to crystalline order. 
In both cases, little information can be extracted about 
the form of the interactions. 

We have derived an expression for the pair correlation 
function in the limit of point-like disorder and low flux 
line density, Eq. (|l5|) . Our theory shows that the small r 
behavior of the pair correlation function depends both on 
the pair interaction and also on the availability of strong 
pinning sites. We found that the distribution of the latter 
is not necessarily the same as the distribution of a pri- 
ori available pinning sites. Rather, the dynamics favor 
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strong pins. Inspection of the stable configurations re- 
vealed that the distribution of dynamically selected pins 
is consistent with a process of choosing the strongest out 
of n available pins. For the parameters used in our sim- 
ulations, n = 2 described the dynamically selected pin- 
ning strength distributions very well. We expect how- 
ever that the parameter n, while apparently independent 
of the width of the pinning strength distribution, should 
depend on both the concentration of pins and also the 
penetration depth of the flux line interaction. A more 
detailed investigation will be left for future work. 

We now turn to a comparison of the qualitative fea- 
tures of the experimentally obtained pair correlation 
functions with those obtained from simulation. We find 
that the width and location of the first peak of g(r) is dif- 
ferent: while the first peak of the pair correlation function 
as obtained from experiments (c/. Fig. is broad and oc- 
curs around the mean flux line spacing, we see a much 
narrower peak in our simulations (cf. Fig. |J). Moreover 
the first peak occurs at a distance significantly smaller 
than the mean flux line spacing. Thus the numerical 
simulations have larger fluctuations in the local flux line 
density than are seen in experiment. A reason for this is 
that our particle-like treatment of flux lines ignores their 
magnetic properties. One expects that including bound- 
ary conditions for a uniform applied external magnetic 
field would favor a homogeneous flux line concentration 
and thus would penalize the large fluctuations in our sim- 
ulations. 
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